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^^ Abstract 

o 

CN This paper analyzes the outage performance in finite wireless networks. Unlike most prior works, which either 

pi I assumed a specific network shape or considered a special location of the reference receiver, we propose two general 

^^ frameworks for analytically computing the outage probability at any arbitrary location of an arbitrarily-shaped finite 

m 

C^ wireless network: (i) a moment generating function-based framework which is based on the numerical inversion 

^ " ^ of the Laplace transforms of probability distribution and (ii) a reference link power gain-based framework which 

r^ 

*; exploits the distribution of the power gain between the reference transmitter and receiver The boundary effects 

I ^1 are accurately accounted for using the exact probability distribution function of the distance of a random node 

,__! from the reference receiver in finite networks. For the case of the node locations modeled by a uniform Binomial 

> 

Aq point process and independent and identically distributed Nakagami-m fading channels, we demonstrate the use of 



o 
en 



X 



the proposed frameworks to analytically evaluate the outage probability at any location in the network of either 
a disk or polygon shape. Our frameworks generalize and reproduce many existing results in the literature. The 
analysis illustrates the location-dependent performance in finite wireless networks and highlights the importance 
of accurately modeling the boundary effects. 
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I. Introduction 

A. Motivation 

Outage probability is an important performance metric for wireless networks operating over fading 
channels [1|. It is commonly defined as the probability that the signal-to-interference-plus-noise ratio 
(SINR) drops below a given threshold. The analysis of the outage probability and interference in wireless 
networks has received much attention recently [|2|-[|9|. For the sake of analytical convenience and tractabil- 
ity, all the aforementioned studies and many references therein assumed infinitely large wireless networks 
and often used a homogeneous Poisson point process (PPP) as the underlying model for the spatial 
node distribution. A homogeneous PPP has a nice property of being stationary i.e., the node distribution 
is invariant under translation. This gives rise to location-independent performance, i.e., statistically the 
network characteristics (such as outage probability) as seen from a node's perspective are the same for all 
nodes. Mathematical tools from stochastic geometry [[T|, [|7|, [[8| have been applied to obtain analytical 
expressions for the outage probability in infinitely large wireless networks. The outage analysis in infinite 
wireless networks has also been extended to wireless networks with stationary but non-homogeneous node 
deployments [jSJ, [|7|, pO| , as well as to coexisting networks sharing the same frequency spectrum jllf- 

In practice, many real-world wireless networks comprise of a finite number of nodes distributed at 
random inside a given finite region. The boundary effect of finite networks gives rise to non- stationary 
location-dependent performance, i.e., the nodes located close to the physical boundaries of the wireless 
network experience different network characteristics (such as outage probability) as compared to the 
nodes located near the center of the network. As a result, the modeling and performance analysis of finite 
wireless networks requires different approaches as opposed to infinite wireless networks. For example, 
when a finite number of nodes are independently and uniformly distributed (i.u.d.) inside a finite network, 
a Binomial point process (BPP), rather than a PPP, provides an accurate model for the spatial node 



distribution [ 15 1, [ 16|. Unlike infinite wireless networks, deriving general results on the outage probability 
in finite wireless networks is a very difficult task because the outage performance depends strongly on the 
shape of the network region as well as the location of the reference receiver. In this work, we would like 
to investigate whether there exist general frameworks that provide easy-to-follow procedures to derive the 
outage probability at an arbitrary location in an arbitrarily-shaped finite wireless network. 



B. Related work 

Since it is difficult to derive general results on the outage probability in arbitrarily-shaped finite wireless 
networks, most prior works focused on a specific shape (such as disk) and computed the outage probability 
at a specific location (such as the center of the network region). A few recent studies presented outage 
characterizations at the center of a BPP network with more general shapes p3| , [ 17|, 1 18 1. Specifically, the 



analytical expression for the moment generating function of the aggregate interference seen at the center 
of a BPP network was presented in [17]. The results were extended for the case of spatial multiplexing 



with the maximal ratio combining and zero forcing schemes in [18|. Using the probability distribution 
function (PDF) of the nearest neighbor in a BPP, a lower bound on the outage probability at the center 
of the wireless network for a simple path loss model was computed in [15|. The exact closed- form 
outage probability in a class of networks with isotropic node distribution (i.e., the node distribution is 
invariant under rotation) was derived in (19). For finite networks, the results in (T9\ can only be applied 
to very limited cases preserving the isotropic node distribution, such as a disk-shaped network with BPP 
node distribution. The outage probability at an arbitrary location of an arbitrarily- shaped finite wireless 



network was studied in [20|. This work focused on deriving closed-form expressions for the conditional 



outage probability, which is conditioned on the locations of all the interfering nodes in the network. For 



the (unconditional) outage probability averaged over the spatial distribution of nodes, the authors in [20| 
presented the analytical result for the special case of an annular- shaped network with the reference receiver 
at the center. For other shapes and receiver locations, the authors in [ |20[ suggested the use of Monte Carlo 
simulations to compute the outage probability. Therefore, it is still largely an open research problem to 
find general frameworks for deriving the outage probability at an arbitrary location of a finite wireless 
network with an arbitrary shape. 

C. Contributions 

In this paper, we focus on the outage probability analysis for a reference transmitter-receiver link in 
the presence of M interferers and additive white Gaussian noise in a finite wireless network. We present 
analytical frameworks for computing the outage probability at an arbitrary location in an arbitrarily-shaped 
finite wireless network. Specifically, we make the following contributions in the paper: 

• We propose two general frameworks for the exact calculation of the outage probability in arbitrarily- 
shaped finite wireless networks in which the reference receiver can be located anywhere. 



- The first framework, named the moment generating function-based (MGF-based) framework, is 
based on the numerical inversion of the Laplace transforms of the probability distributions. It is 



inspired from the mathematical techniques developed in pT[ |, p2| and is valid for any spatial 
node distribution and any fading channel distribution. To the best of our knowledge, this is the 
first time that such an approach has been applied in the context of finite wireless networks. 
- The second framework, named the reference link power gain-based (RLPG-based) framework, 
exploits the distribution of the power gain between the reference transmitter and receiver. It is 



based on the combination and generalization of the frameworks proposed in |17|, pO| and is 



valid for any spatial node distribution and a general class of fading channel distribution proposed 



in [231 



In order to demonstrate the use of the proposed frameworks, we consider the case where the nodes are 
independently and uniformly distributed in the finite wireless network, with node locations modeled 
by a uniform BPP. The fading channels between all links are assumed to be independently and 
identically distributed (i.i.d.) according to a Nakagami-m distribution. We use the exact probability 
distribution function of the distance of a random node from the reference receiver in uniform BPP 
networks [ |16[ to accurately capture the boundary effects in the two frameworks. 

- We demonstrate the use of the two frameworks in evaluating the outage probability for the 
important case of a reference receiver located anywhere in a disk region. We show that the 
known outage probability results in the literature for the disk region arise as special cases in our 
two frameworks. 

- We further consider the case of an arbitrary located reference receiver in a convex polygon and 
present an algorithm for accurately computing the outage probability. This fundamentally extends 
the prior work dealing with the outage probability analysis in finite wireless networks. 

- We show that due to the boundary effects, the outage probability in regular polygon regions 
does not provide any meaningful bounds for the outage probability in an arbitrarily-shaped finite 
region. This highlights the importance of the proposed frameworks, which allow accurate outage 
probability computation for arbitrarily- shaped finite wireless networks. 



D. Notation and Paper Organization 

The following notation is used in this paper. Pr(-) indicates the probability measure and E{-} denotes 
the expectation operator, i is the imaginary operator and Re{-} denotes the real part of a complex- valued 
number. |^| denotes the area of an arbitrarily-shaped finite wireless network A. V{x) = J^ t^~^ exp(— t)(it 
and r(a,x) = J^ t^~^ exp{—t)dt are the complete gamma function and incomplete upper gamma func- 



tions, respectively [24|. 2-^"! (■, ■, ■, ■) is the Gaussian or ordinary hypergeometric function [24|. In the 
context of random variables (RVs), upper case italic letters (e.g., Z) represent random variables while 
lower case italic letters (e.g., z) represent deterministic variables or constants, fzi^) and Fz{z) denotes 
the probability distribution function and the cumulative distribution function (CDF) of a random variable 
Z. Cz{s) = E{exp(— sZ)} indicates the Laplace transform with respect to the random variable Z. Also, 
Table |l] summarizes the main mathematical symbols used in this paper. 

The remainder of the paper is organised as follows: The system model and assumptions are presented 



in Section |IIj The proposed general frameworks are described in detail in Section III For the case of the 



node locations modeled by a Binomial point process and i.i.d. Nakagami-m fading channels, the evaluation 



of the proposed frameworks for arbitrarily- shaped finite wireless networks is detailed in IV The outage 



probability analysis for the case of a reference receiver located anywhere in a disk and a convex polygon 



is presented in Section M and Section VI respectively. The derived results are used to study the outage 



probability in Section VII Finally, the paper is summarized in Section VIII 



II. Problem Formulation and System Model Assumptions 

Consider a wireless network with M + 2 nodes which are located inside an arbitrarily- shaped finite 
region ^ G M^, where M^ denotes the two-dimensional Euclidean domain. The M + 2 nodes consist of a 
reference transmitter Xq, a reference receiver Yq and M interfering nodes. The M interfering nodes are 
distributed at random within the region A. Throughout the paper, we refer to Xj (z = 1, 2, . . . M) as both 
the random location as well as the zth interfering node itself. The reference receiver Yq is not restricted 
to be located at the center of the finite region but can be located anywhere inside the region A. The 
reference transmitter Xq is assumed to be placed at a given distance tq from Yq. Let i?j (z = 1, 2, . . . M) 
denotes the Euclidean distance between the ith interferer Xj and the reference receiver Yq. 

We focus on the performance of the reference link comprising the reference transmitter Xq and the 
reference receiver Yq, in the presence of M interfering nodes and noise. We consider a path-loss plus 
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Number of interfering nodes 

Arbitrarily-shaped finite region 

Radius of disk or circum-radius of a polygon 

Reference transmitter and its location 

Reference receiver and its location 

Euclidean distance between Xo and Yo 

Transmit power for reference transmitter Xo 

Transmit power for interferer Xi 

Path-loss coefficient 

Nakagami-m fading parameter for reference link 

Nakagami-m fading parameter for interference links 

Additive white Gaussian noise (AWGN) power 

Signal-to-noise ratio 

SINR threshold 

Outage probability 
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ith interfering node and its location 
Euclidean distance between Xi and Yo 
Power gain for the reference link 
Power gain for the ith interference link 
Aggregate interference 
Signal-to-interference-plus-noise ratio (SINR) 



block-fading channel model. Let Go represent the instantaneous power gain for the reference link and 
Gi represent the instantaneous power gain between X,; and Yq. We assume all interferers Xj are in the 



far-field of the reference receiver Yq to avoid the effect of the singularity in the path-loss model [25|. 
Thus, the path-loss function can be expressed as 



/(r) 



(1) 



where r denotes the distance and a is the path-loss exponent, which typically lies in the range 2 < a < 



6 [|26|. 
Let Po and P, denote the transmit powers for Xo and Xj, respectively. The aggregate interference power 



at the reference receiver Yq is then given as 

M 
I = Y,P^G^Ri". (2) 

The instantaneous signal-to-interference-plus-noise ratio 7 at the reference receiver Yq is given by 

_ PoGoTq " _ Gq 

^ ~ N + I ~ i + — ?— ' ^^ 

po ^ Pon7" 

where po = (-Po^(7°)/^ is defined as the average signal-to-noise ratio (SNR) and A^ is the additive white 
Gaussian noise (AWGN) power. 

We characterize the performance of the reference link in the presence of the aggregate interference and 
AWGN by using the outage probability. An outage is said to occur when the SINR 7 falls below a given 
SINK threshold (3, i.e., 

e = Pr(7</3). (4) 

We are interested in obtaining the average outage probability in the arbitrarily- shaped finite wireless 
network after un-conditioning with respect to the spatial node distribution and the fading distribution. 
This is addressed in the next section. 

III. Proposed Frameworks 

In this section, we propose two analytical frameworks to compute the outage probability in arbitrarily- 
shaped finite wireless network. We also discuss the merits and demerits of each framework. The first 



framework, named the moment generating function-based (MGF-based) framework, is inspired from [21 1, 



. The basic principle of this framework is the accurate numerical inversion of the Laplace transform 



of probability distribution for an appropriately defined random variable [21 1. The second framework. 



named the reference link power gain-based (RLPG-based) framework, is based on the combination and 



generalization of the frameworks proposed in [17|, [20|. The basic principle of this second framework is 



to find the cumulative distribution function of the reference link's power gain, which can then be used to 
find the outage probability. These frameworks are discussed in detail in the following subsections. 

A. Moment Generating Function-based (MGF-based) Framework 

In this framework, it is necessary to define a suitable random variable. Substituting (|3]) into (|4]) and 
rearranging, we have 



We define a random variable Z as 



Hence, (|5]) can be re-written as 






e = Pr(Z>r') = l-i^z(r'), (7) 



where Pr(-) denotes the probability measure and -Fz(-) is the CDF of the random variable Z. 

In general, it is not possible to obtain a closed-form solution for Fz{/3^^). Hence, we use numerical 



inversion of Laplace transform to find Fz{l3 ). Following [21 1, we have 



^ 6=0 ^ / c=0 "^ *- ■' 

where D^ = 2 (if c = 1) and D^ = 1 (if c = 2, 3, . . .), s = {A + i27vc)/{2/3-^) and Re{-} denotes the real 
part. The three parameters A, B and C control the estimation error. The selection of the values for A, B 
and C for accurate numerical inversion will be discussed later in Section IVII-Al 

Using the definition of the Laplace transform of a random variable, we can express Cz{s) as 

where E{-} denotes the expectation operator. 

Combining ([8]) and (|9]), we have the general outage probability expression resulting from the MGF-based 
framework as 






6=0 ^ ^ c=0 




B. Reference Link Power Gain-based (RLPG-based) Framework 

As highlighted earlier, this framework relies on the cumulative distribution function Fcf^igo) for the 
reference link's power gain Gq. We adopt the following general modeQfor the CDF [ 23 1 



Fcoigo) = 1 - X] exp{-ngo)^ankgo, (H) 

neJV k&K 

where the finite sets J\f, /C C N (where N is the set of natural number) and a„fc are the coefficients. 
With the proper choice of a„jt, J\f and /C, Fcoigo) can represent different types of distributions for the 



power gain [23 1. For example, when Af = {1}, /C = {0} and a^fc = 1, then Fcoigo) = 1 — exp(— (70) 



We note that it may be possible to consider other general classes of fading channels 27 - 29 , but this is outside the scope of this paper 



and consequently /Go(fl'o) reduces to the Exponential distribution, which corresponds to the reference 
link undergoing Rayleigh fading [30|. Furthermore, if A/" = {mo}, /C = {0,...,mo — 1} and ank = 



f, then Fgo{9o) = {9o^° "^cT" 6xp(— mo5'o))/r(mo) and consequently faoido) reduces to the Gamma 



distribution [31 1 which corresponds to the reference link experiencing Nakagami-m fading. 



We proceed by re-writing the outage probability in (|4]) as 

Go 



e = E/ <^ Pr 
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Ej{Pi{Go<f3 



1 



Po Poro" 



(12) 



po Poro" 

where E/{-} denotes the expectation with respect to the aggregate interference. Using Fcoigo) shown in 



(11), (fT2b can be written as 



e = E, <^ F, 



'Go 



/3 , /3 
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1 - V exp ( -n— ) V ank'Ei < exp ( -n-—^^I I I — + „ _ 
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(13) 



where the aggregate interference / depends on the power gain distribution and the distance distribution. 
Hence, ( [T3| ) can be further expanded in terms of G, and Ri. First, we focus on the expansion of the term 
^/j f — + j^-^Ij . Based on the binomial theorem [32| and using (|2]), we have 



exp I —n 



exp —n 



/3 



/3 



/3 
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k-j 



/3 



Poro 



M 
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(14) 



,i=l 



^A/ 



Then, following the multinomial theorem [33 1, we can further expand the term ( X]j=i PiGiRi " ) into 

^ M 



,i=l 



tl+t2 + ...+tM=J ^ ' ' ' / j=l 



A/ 



(15) 



where U {i = 1,2, ...M) is a non-negative integer and the multinomial coefficient (^ t '^ t ) ^ t \t ! ' t ,\ - 
Combining ([14]) and ( [T5| ) and substituting back into ( [T3] ), we have the general outage probability 
expression resulting from the RLPG-based framework as 



e =1 - E exp (-n- j E«"^E ( , 



E 



G„R, < exp 



-72- 



>oro- 



M 



Po 

Af 



fc-j 



/3 



^oro-' 



E 



tl+t2---+tA/=i 



J 
tl, ^2, ..., t 



A^ 



X 



-E^Ai?r n(^^^^^^"")' 



(16) 



i=\ 



i=\ 



C. Comparison of the Two Frameworks 



The proposed frameworks have their own merits and demerits. The two general formulations in (10) 



and ( [T6| ) are valid for any spatial node distribution in an arbitrarily-shaped finite wireless network and any 
location of the reference receiver inside the arbitrarily- shaped finite region. The MGF-based framework 
expression in pO]) is valid for any fading channel distribution. The RLPG-based framework expression 



in ([16]) is valid for a general class of fading channel distribution expressed in ( [TT] ). There are further trade- 



offs in the computational complexity of (10) and (16) and their ability to yield closed- form analytical 



expressions in certain cases of interest. This will be further discussed in Remark [T] in Section IV and in 
Section Ivn^ 



The evaluation of the outage probability in ([10]) and ( [T6] ) requires the knowledge of the joint probability 



distribution function of the distance Ri and the power gain Gi. In the next section, in order to demonstrate 
the evaluation and use of the proposed frameworks, we consider the case of that the nodes are distributed 
at random inside an arbitrarily- shaped finite wireless network according to a BPP and all fading channels 
are i.i.d. Nakagami-m fading channels. 

IV. Outage probability for BPP and NAKACAMi-m Fading Channels 
In this section, we consider the scenario where: 
Al. The M interfering nodes are independently and identically (i.i.d.) distributed at random inside an 
arbitrarily- shaped finite wireless network, i.e. the nodes are distributed at random according to a 
general BPP yj Definition 2.12]. 
A2. The fading channels are independently and identically distributed (i.i.d.). 
A3. The transmit powers Pq and Pj are normalized to unity. 

A4. The fading channels follow a Nakagami-m distribution. The Nakagami-m distribution is widely used 
in the literature to model the distribution of the signal envelopes in various fading environments. 



such as the land-mobile and indoor-mobile multipath propagation environments [26|. The parameter 
m, which lies in the range 1/2 to oo, describes the severity of the fading channel. Note that it 
is not necessary for m to be an integer number only. When m < 1, the Nakagami-m distribution 
provides a close approximation to the Nakagami-g (Hoyt) distribution with parameter mapping m = 
((l+g^)^)/(2(H-2g^)). Additionally, when m > 1, the Nakagami-m distribution closely approximates 
the Nakagami-n (Rice) distribution with parameter mapping m = ((1 + n^)^)/(l + 2n^) []26|. It is 
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well known that m = oo corresponds to the no-fading case, m = 1 represents the special case of 



Rayleigh fading and m = 1/2 represents that unilateral Gauss distribution [26 1, which corresponds 

to the most severe Nakagami-m fading. 
A5. The nodes are distributed at random according to a uniform BPP [[I] Definition 2.11]. This means 

that the nodes are independently and uniformly distributed (i.u.d.) inside the arbitrarily-shaped finite 

wireless network. 
As a consequence of assumptions Al and A2, the joint PDF of the distance Ri and the power gain 
Gi can be decomposed into the individual PDFs, which are denoted as fRX^i) and fdigi), respectively. 
Due to assumption Al, the distribution of Ri is the same for all i. Similarly, due to assumption A2, the 
distribution of Gi is the same for all i as well. Thus, we can drop the index i in Ri, Gi, fuX'^i) and 
fcM and let //{.(n) = fnir) and fcXQi) = fcig)- 

Let fcoigo) denotes the PDF of power gain for the reference link. From assumption A4, since the 
fading coefficients for both the reference link and the interference links are modeled using a Nakagami-m 
distribution, i.e., fcoigo) = fcig), the distribution for the power gains Go and G can be modeled as a 



Gamma distribution with the following PDF [31 1 



„m-l,^m 

fcig) = -p. X exp(-m^). (17) 

1 [m) 

Note that we will represent the Nakagami-m fading parameter for the reference and interfering links as 

mo and m, respectively. Even though the power gain distributions are identical, we will still use mo and 

m to distinguish the reference link from the interference link for the sake of analytical convenience. 

A. MGF-based Framework 

Using assumptions Al— A3, i.e., i.i.d. random nodes and i.i.d. fading channels, the general outage 



probability expression in ([TO]) for the MGF-based framework can be simplified to 



2-^ exp(l) J^ fB\ '^ (-ir I ^^'^ i'"P l"^i l^^'^ i'"P l^" 
1^\ hi 1^ n ^^ 



' fe=o ^ ^ c=0 "^ 

(18) 
where EgoI-} represents the expectation with respect to Gq and 'Ejc,r{-} represents the expectation with 
respect to R and G. 
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Using assumption A4 and substituting ( [17] ) in ( [T8| ), the outage probability in ( fTS] ) can be expressed as 

2-^exp(f) 



e=l- 



Re 



■A\ B /^x C+b ,_^y 



13- 



:?^E : E 



b=0 ^ ^ c=0 ^ 



PCX) 

j exp 



Po5oy v^^'^i^'^Pv r-rco 



■X 



-sGiJ- 



r(mo) 



^^2 ^^^bJI exp{-mogo)dgo 



where the expectation in ( [T9| ) can be expressed as 
Eg,r < exp 



To'^Go 



m'" I m + ^ 



fR{r)dr. 



Note that mo and m in (19) can take any values (whether integer or non-integer). 



>, 



(19) 



(20) 



B. RLPG-based Framework 

Using assumptions Al— A3, i.e., i.i.d. random nodes and i.i.d. fading channels, the general outage 



probability expression in ( [T6| ) for the RLPG-based framework can be simplified to 

k /, \ / n \ k-j ^ . s M 

e=i-> expi -n- ) 7 a„fc7 ^( . ) ( - J ^pr^y ^ [. .^ . ]\{^g,r{^u} 



1 - E -p ( -«£) E - 1 G) O" ""»''■ - v«. *. .... <M. 

^'^^ fce/c i=o ^"'^ ^^'^^ ti+t2...+tM=j ^ ^' "^' ' ^^-'z j=i 



n&Af 



(21) 



where fit- = exp {—mQ(3r'^GR °) (Gi? ") * is defined for analytical convenience. 



Using assumption A4, the outage probability in ( [2T| ) can be expressed for the case of Nakagami-m 
fading with integer niQ as 



.=.-.-.3£-|i:(:)0'-«.'...?.,.,(i: 



EcRJfitJ, (22) 



where the expectation in ( [22j ) can be expressed using ( [17] ) as 



EG,i?{fitJ 



mi'^-Q^i l'^ _^ /:J^a^^^-Q^I ™"*^r(^ + t. 



R 



Tim) 



-fR{r)dr. 



(23) 



Remark 1: Both frameworks have their merits and demerits. For the MGF-based framework, the outage 



probability in ( 19) involves a double integration (an integration with an expectation term in the integrand). 



In general, it is not possible to obtain a closed-form for the integration part in ( ]T9] ) because the expectation 
term is raised to a power factor of M (M > 2). However, in certain cases, the expectation term can be 



expressed in closed-form. For the RLPG-based framework, the outage probability in ( [22] ) involves a single 



integration which admits closed-form results in a much larger number of cases. However, m^ in ([22]) is 



constrained to take integer values only while m can take any value. This is in contrast with (19) where 
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both mo and m can take any (integer or non-integer) value. Thus, we recommend using the RLPG-based 
framework for scenarios with integer niQ. However, there are scenarios, i.e., non-integer mo, where only 
the MGF-based framework can be applied. 



C. Summary 



Summarizing, ( |T0| ) and ( |T6| ) take the form of ( |T9| ) and ( [22| ), respectively, for M interfering nodes i.i.d. at 
random inside an arbitrarily- shaped finite wireless network with i.i.d. Nakagami-m fading channels. 



Both ( [191 ) ^iid ([22]) require the knowledge of the distance distribution fnir), i.e. the PDF of the 
distance of a random node from the reference receiver Yq. The distance distribution //j(r) is dependent 
on the underlying random model for the node locations. For the uniform BPP (assumption A5), which is 



considered in this work, the distance distribution //j(r) is derived in [15| for the special case when the 
reference receiver is located at the center of a convex regular polygon. Recently, the result in [15] was 



generalized in [16| for the case when the reference receiver is located anywhere inside a convex regular 



polygon. It was shown in [ 16| that for an arbitrary location of the reference receiver inside convex regular 
polygon, the distance distribution fR{r) can be a complicated piece-wise function because of the boundary 
effects. We note that the approach in [[T6| is also applicable for arbitrarily-shaped convex polygons (see 
Appendix |A] for details). Once /^(r) is given, both ( [T9| ) and ([22]) and can be accurately evaluated. 



In the next two sections, we show how ( [19] ) and ( ]22[ ) can be evaluated in disk and polygon regions, 
which are commonly used in the literature for the modeling of wireless networks. We also illustrate how 
the proposed frameworks can be applied in the case of arbitrarily-shaped finite wireless networks with 
arbitrary location of the reference receiver. 



V. Outage Probability in a Disk Region 



Consider the scenario that the region ^ is a disk of radius IZ, as shown in Fig. [T] The reference receiver 
is assumed to be located at a distance d from the center of the disk. Then, the distance distribution fR{r) 



can be exactly expressed as [T6| 

fR{r) = 



2r_ 

-^ arccos 



0<r <n-d: 



{"^^tr^)^ n-d<r<n + d- 



(24) 
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Fig. 1. Illustration of a finite wireless network with arbitrary location of reference receiver Yo in a disk region of radius ??.(+ = center 
of disk, • — interfering node, A = reference receiver, ■ — reference transmitter). 



Note that substituting d = in ( |24l ) gives the distance distribution for the special case that Yq is located at 



the center of the disk. Similarly, substituting d = TZ in (|24]) gives the distance distribution for the special 
case that Yq is located anywhere on the circumference of the disk. 



A. MGF-based Framework 



Substituting (24) in ([20]), we find that the expectation has a closed-form only for the first part of the 



range (0 < r < 7^ — rf) in ( [24] ). For the second-part of the range {IZ — d < r <1Z + d), the integration 
does not have a closed-form due to the arccos(-) function. The result is shown in ( [25] ). 



Ecij < exp 



-sGR- 



rn-°a 



'-^o 



2m^g^ {JZ - d) 
"^^(2Tam)7rp) 



2+am 



iFi 



m, 



+ m, 1 



TTl, 



a 



a 



gprn {TZ - df 



+ 



2mJ 



TT 



7^2 



n+d 



Tl-c 



m + 



9o 



V2 + rf2_^2~ 

r arccos | ; 1 dr. 

2dr 



(25) 



Substituting ([25]) in ( 19 ), the outage probability can be numerically evaluated using a standard mathematical 



software package such as Mathematica. This will be discussed in detail in Section VII-A 



B. RLPG-based Framework 



Substituting ( [24] ) in ( [23] ) and after some manipulations, we get 

2 



■CJG,fl|"tiJ = ^o/^ , r;^7 — ^ 2-ri 



+ 



7^2(2 + am)T{m 



TiTVV{m) U 



2 m(n - dV 

- + m,m + ti,H \-m, -— 

a a pr^mo 

/ 2 I j2_-i-,2~ 

(r ) ' (^m + pTq uiQr j r arccos I :r-. I dr. 



2dr 



(26) 
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In general, the integration in ([26]) also does not have a closed-form due to the arccos(-) function. It is 



possible to use the Gauss-Chebyshev integration technique [34| in order to obtain an approximate closed- 
form expression. However, in our investigations, we found that a summation over a large number of terms 
(> 1000) was required in our case for accurate evaluation. Hence, we do not pursue approximations and 



instead directly substitute ([26]) in ( [22] ) to obtain the outage probability. 

Special case: For the case of the reference receiver Yq located at the center of the disk, d = 0. Substituting 



this value in ([26]) and then substituting the result in ( 122] ), the final expression for the outage probability 
simplifies to 



exp 



^'J/ fc=0 j=o ^-^ 



M 

n 

i=l 



2m™' (/3rQ mo 



,—m~ti 



7^"'"^('m + ti) 



(2 + am)T{m) 



Po 



iFi 



(/3r, 



o) 



E 

tl+t2...+tni=j 

2 2 

— \-m,m + ti,l-\ — 
a a 



J 

ti,t2, ...,tM 

m, 



X 



(3r^mo 



(27) 



C. Comparison with Existing Results for a Disk Region 

The outage probability for a finite number of nodes i.u.d. in a disk region has been widely considered in 
recent literature and we show how existing results can be obtained as special cases from our frameworks. 



For the MGF-based framework, with mo = 1, results from ( [19] ) and ( [25] ) is equivalent to the result 



in p7] eq.(24)]. For the RLPG-based framework, with the reference receiver located at the center of the 



network, (127]) is identical to the results in 120 eq.(44)]. Finally, the closed-form outage probability in 



an isotropic disk region is presented in 119|. Their results, which are limited to the special case of i.i.d. 
Rayleigh fading channels with path loss coefficient a = 2 and 4 only, match exactly with the results 
from ( l22] ) and ( |26] ). This will be further demonstrated in Section fVIIj 



VI. Outage Probability in Polygons and Arbitrarily-Shaped Regions 

In this section, we illustrate the exact computation of the outage probability in both regular and 
arbitrarily-shaped convex polygon regions. We will consider the following two cases (i) reference receiver 
Yq located at the center of a regular L- sided polygon and (ii) reference receiver Yq located at an arbitrary 
location in an arbitrarily- shaped region. 
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A. Center of Polygon 

Consider the finite region ^ to be a regular L-sided convex polygon which is inscribed in a circle of 
radius IZ. Then, the area and the interior angle between two adjacent sides are given by 



1^1 



e 



7r(L-2) 



(28a) 
(28b) 



In general, polygon regions are non-isotropic. Hence, there is no single expression for the distance 
distribution //^(r) for an arbitrary location of the reference receiver inside a convex regular polygon. 
For the special case that the reference receiver Yq is located at the center of an L-sided convex regular 



polygon, the distance distribution fR{r) can be expressed as [15|, [16| 

1 I 27rr, < r < 7^sin(f); 



fR{r) 



\A\ 



27rr - 2Lr A, 7^ sin (|) < r < 7^; 



(29) 



where A = arccos ( — z^^ 1 • Using ([29]), we illustrate the computation of the outage probability using 



the two frameworks. 



MGF-based Framework: Substituting (29) in (20) and after some manipulations, we get 



Eg,/? < exp 



-sGR-'' 



27rm'"c^™7^2+^ 



\A\{2 + am){r^s) 
2LmJ^ '■^ 



2^1 



2 2 gomVf 

m, — \- m,l -\ \- m, 



a 



1^1 ^sm(f) 



m + 



9o 



a 



r arccos 



'^sin(f) 



dr. (30) 



Substituting ( pO] ) in ( [T9| ), the outage probability can be evaluated. 

RLPG-based Framework: Substituting (|29]) in ( [23] ), and after some manipulations, we get 

2 2 

- + m, m + tj, 1 H h m, - 

a a 



Eg,_R, |"t,| = 1 t\rc^ . TT^T ^ 2-ri 



\A\{2 + am)V{m) 



mW 



\A\T{m) U 



(r-")*' (m + /3r^mor-") 



-m—ti 



r arccos 



(3r^mo 

nsm{l) 



dr. 
(31) 



Substituting (31) in (22), the outage probability can be evaluated. 



B. Arbitrarily-Shaped Region 

We consider an arbitrarily- shaped convex polygon region as shown in Fig.|2} with side lengths 5*1 = 27^, 
5*2 = \/27l, 3^ = 71 and 5*4 = 7?.. Suppose that Yq is located at the vertex V3. Then, following the derivation 
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Fig. 2. Illustration of an arbitrary location of a reference receiver in an arbitrarily-shaped finite wireless network, with side lengths Si = 27?., 
52 ~ \/2TZ, S'i —TZ and S4 — TZ and vertices Vi, V2, V3 and V4. The areas B2 (shaded in red lines), B3 (shaded in green lines) and C3 
(intersection of red and green lines) are defined in Appendix [aI (• = interfering node, A = reference receiver, ■ = reference transmitter). 



in Appendix \A\ fR^r) can be expressed as 



fnir) 



\A\ 



4 ' ' 



37r„ 



< r < 7^; 

- 3r arccos (f ) , n<r < ^271; 



(32) 



Using ( [32| ), we illustrate the computation of the outage probability using the two frameworks. 

s 
gom{V2ny 



MGF-based Framework: Substituting ([32]) in ( [20] ), the expectation can be expressed as 

2+am 



EcR { exp 



-sGR-' 



37rm™^™(y27^)' 



A\A\{2 + am){r^s) 

'■V2TI 

m 



iFi 



3m" 



90 



2 2 

m, — h m, 1 H \- m, 

a a 



r arccos I — \ dr. 



(33) 



Finally, the outage probability for the case of the reference receiver located at vertex V3 can be evaluated 



by substituting ( ]33] ) in ( ]T9] ). 
Framework: Si 



RPLG-based Framework: Substituting ([32]) in ( 123] ), the expectation can be expressed as 



(V27^)^ 

4|^|(2 + am)r(m 



2 2 

— h m, m + tj, IH h m, 



m 



{V2ny 



fn 



(r ")*' (m + PrQ-mor "^ "* ' r arccos I — \ dr. 



/3r^mo 
(34) 



|^|r(m) Jn ~ ' ' " ' \r 

Finally, the outage probability for the case of the reference receiver located at vertex V3 can be evaluated 



by substituting ^^ in ([22]). 

Summarizing, the procedure for deriving the outage probability for i.u.d. nodes in an arbitrarily-shaped 
convex polygon region with i.i.d. Nakagami-m fading channels is summarized in Algorithm [T] 
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Algorithm 1 Proposed Algorithm 



Step 1: Choose the location of the reference receiver inside the arbitrarily- shaped convex polygon 

region. 

Step 2: Determine fnir) based on the approach summarized in Appendix [A] 

Step 3: Depending on the value of mo, select the appropriate framework to calculate the outage 

probability. 

if rriQ is non-integer then 



Use the MGF-based framework. Substitute /^(r) in ( [20] ) and then ( [20] ) in ([19]) to compute the outage 
probability. 
else 



Use the RLPG-based framework. Substitute fnir) in ( [23] ) and then ( [23] ) in ( |22l ) to compute the outage 
probability. 
end if 

VII. Numerical and Simulation Results 

In this section, we first address the computational aspects of the two frameworks. We also briefly describe 
the Monte-Carlo simulation setup which is used to demonstrate the accuracy of the proposed frameworks. 
We then study the outage probability performance of arbitrarily- shaped finite wireless networks and discuss 
some important insights concerning boundary effects in finite wireless networks. 

A. Computational Aspects of the Frameworks 

We use Mathematica to compute the results for the two frameworks in this paper, with the final plotting 
being done in Matlab. 



For the MGF-based framework, the outage expressions in (10), ([18]) and (19) are a summation over a 



finite number of terms. The number of terms controls the accuracy. Following the guidelines established 



in [ ]22| , we set A = 10 In 10, B = lA, C = 8. This achieves an estimation error of 10^ in most cases. For 
the accurate calculation of the integrations, we use the Nlntegration[ ] command in Mathematica with 
the following settings: PrecisionGoal=10, WorkingPrecision=30, MaxRecursion=20 and Method -> 
ExtrapolatingOscillatory. Further details on these options for the Nlntegration[ ] command are available 
in p5| . It must be noted that the numerical evaluation of double integrations is widely practiced in the 
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TABLE n 

Values of the Main System Parameters. 



Parameter 


Symbol 


Values 


Number of interfering nodes 


M 


10 


Euclidean distance between Xq and Yo 


ro 


5 


Disk radius 


7^ 


100 


SINR threshold 


P 


dB 


Path-loss coefficient 


a 


2, 2.5, 3, 4 


6 


Signal-to-noise ratio 


pa 


20 dB 


Nakagami-m fading parameters 


mo = m 


0.5, 1, 1.5, 


2 



wireless communications literature p6| |. 

For the RLPG-based framework, if both itiq and M are large, then the computation of all possible 



integer results for tj in ti + ^2 + ••• + ^m = j (j = 0, ...,mo — 1) in ( [76] ), pT] ) and ( [22] ) can be time- 
consuming. This is due to the fact that we need to use M for loops to find the complete results. However, 
when either mo or M is a small number we can pre-compute these results, as suggested in [|20|, and store 
them as a matrix for use in computations. 



B. Parameter Settings and Simulation Setup 

We use Monte Carlo simulations in Matlab to generate the simulation results. For the simulation results. 



we have used the acceptance-rejection method [37| to uniformly distribute the users inside an arbitrarily- 
shaped region and averaged the results over 1 million simulation runs. Unless specified otherwise, the 
values of the main system parameters are summarized in Table [Ilj 



C. Disk Region 

Fig. [3] plots the outage probability, e, versus the distance of the reference receiver from the center 
of disk, d, for path-loss coefficients a = 2,3,4,6 and i.i.d. Rayleigh fading channels. The solid lines 



are plotted using ([22]) and (26), i.e., the RLPG-based framework. We compare these results with (i) the 



results from the MGF-based framework, obtained using ( [T9| ) and ( [25] ), (ii) the results from p9| which are 
applicable for a = 2,4 only and (iii) Monte Carlo simulation results. For the sake of clarity, the Monte 
Carlo simulation results for shown for a = 3, 6 only. We can see that the simulation results match perfectly 
with the our analytical results, which is to be expected since we are evaluating the outage probability 
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Fig. 3. Outage probability, e, versus tiie distance of tiie reference receiver from the center of disk, d, for Ad = 10 interferers i.u.d. in a disk 
of radius 7?. — 100, with reference link distance ro = 5, path-loss coefficients a = 2, 3, 4, 6, i.i.d. Rayleigh fading channels {mo = m = 1), 
SINR threshold /? = dB and SNR po = 20 dB. 



exactly. In addition, the results in [19| also match perfectly with our results. These comparisons verify 
the accuracy of the proposed frameworks. We can see from Fig. |3] that, for all values of the path loss 
coefficient, the minimum value of the outage probability occurs when the reference receiver is located at 
the circumference. Additionally, the outage probability is almost constant as long as the distance of the 
reference receiver from the center of disk is approximately less than IZ — tq. This helps to illustrate the 
boundary effects in finite wireless networks. 

Fig.|4]plots the outage probability, e, versus the number of interfering nodes, Af , for path-loss coefficients 
a = 2,3,4,6 and i.i.d. Rayleigh fading channels. The results are plotted using pTj). As expected, the 



outage probability increases as the number of interferers increases. In addition, as the path loss coefficient 
increases, the aggregate interference power is reduced and thus the outage probability decreases, e.g., for 
a target outage probability e = 10^^, the network can tolerate only 7 interfering nodes for a = 2 , which 
increases to 32 nodes for a = 6. 

Fig. [5] plots the outage probability, e, versus the SNR, p^, for i.i.d. Nakagami-m fading channels 
(mo = rn = 0.5, 1,1.5, 2), path-loss coefficient a = 2.5 and the reference receiver located at the center 
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Fig. 4. Outage probability, e, versus the number of interfering nodes, M, i.u.d. in a disk of radius 7?. = 100, witli reference link distance 
ro = 5 and reference receiver located at the center of the network, path-loss coefficients a = 2, 3, 4, 6, i.i.d. Rayleigh fading channels 
(mo = m = 1), SINK threshold ^ = dB and SNR po = 20 dB. 



of the disk. The solid lines are generated using the MGF-based framework ((|T9|) and ([25])), since in this 
case we have mo = m taking non-integer values. For rriQ = m = 1.5 we also plot the arithmetic and 
the geometric means using the outage probabilities for niQ = m = 1 and mo = m = 2, respectively. We 
can see that the arithmetic and geometric means do not match the exact value of the outage probability. 
In addition, the result for niQ = m = 1 does not provide a tight bound on the outage probability when 
mo = m < 1. These issues highlight the importance of the MGF-based framework since the RLPG-based 
framework cannot handle non-integer mo values. 



D. Polygon Region 

Fig. [6] plots the outage probability, e, versus the SNR, po, for equilateral triangle, square, hexagon and 
disk shaped regions inscribed in a circle of radius TZ = 100, reference receiver located at the center of the 
network, path-loss coefficient a = 2.5, i.i.d. Nakagami-m fading channels with mo = m = 3. The results 
are plotted using the RLPG-based framework (([22]) and ([3T])). We can see that the outage probability is 
lowest for the disk region and highest for the equilateral triangle. This is because for the same circum- 
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Fig. 5. Outage probability, e, versus the signal-to-noise ratio, po, for i.i.d. Nakagami-m fading channels and mo = m = 0.5, 1, 1.5, 2, with 
M = 10 interferers i.u.d. in a disk of radius 7?. — 100, reference link distance ro = 5 and reference receiver located at the center of the 
network, path-loss coefficient a — 2.5 and SINR threshold 13 — dB. 



radius, the equilateral triangle has the smallest area. Therefore, the density of the interferers increases and 
hence the outage probability also increases. This emphasizes the importance of accurately characterizing 
the finite area networks. 

In order to illustrate and investigate the boundary effects in finite wireless networks, we plot the 
ratio of the outage probability for the reference receiver located at the center and the corner of an 
L = 3,4, 5, 6, 7, 8, 9-sided polygon for the case of no fading (mo = m = oo) and i.i.d. Rayleigh fading 
(mo = m = 1), with path-loss coefficient a = 2.5. The results are plotted using the MGF-based framework. 
For the case of the reference receiver located at the center of a polygon, we use the distance distribution 
fnir) in ([29]). For the case of the reference receiver located at the corner of a polygon, the distance 
distributions /ij(r) can be determined using the procedure summarized in Appendix |a] and are not included 
here for brevity. The figure shows that for an equilateral triangle (L = 3), the ratio is the highest, which 
represents the most severe boundary effect. As L increases and the shape of the polygon approaches a disk, 
the boundary effect is relatively less severe. Comparing the curves for mo = m = oo and mo = m = 1, 
for the low node density, we can see that i.i.d. Rayleigh fading helps to reduce the boundary effects on 
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Fig. 6. Outage probability, e, versus the signal-to-noise ratio, po, for equilateral triangle, square, hexagon and disk shaped regions inscribed 
in a circle of radius TZ = 100, with M = 10 interferers i.u.d. in each region, reference link distance ro = 5 and reference receiver located 
at the center of the network, path-loss coefficient a = 2.5, i.i.d. Nakagami-m fading channels and mo — m — 3 and SINR threshold /3 = 
dB. 

the outage probability to some extent (the curve for uiq = m = 1 is lower in this case). However, the 
boundary effect is still quite severe for popular finite network shapes such as square and hexagon. 



E. Arbitrarily-Shaped Finite Region 

Fig. [8] plots the outage probability, e, versus the SNR, po, with arbitrary location of reference receiver 
in the arbitrarily- shaped finite region in Fig. 2, i.i.d. Nakagami-m fading channels with uiq = m = 3 
and path-loss coefficient a = 2.5. We consider and compare the following scenarios for the reference 
receiver located at: (i) vertex V3 in the arbitrarily- shaped finite region shown in Fig. 2, (ii) vertex V2 in 
the arbitrarily- shaped finite region shown in Fig. 2, (iii) vertex of an equilateral triangle, (iv) vertex of 
a square and (v) circumference of a disk. The equilateral triangle, square, disk regions are assumed to 
have the same area as the area of the arbitrarily- shaped finite region in Fig. 2, i.e., |^| = 15000. The 
results are plotted using the RLPG-based framework, using the distance distribution results summarized 
in Appendix |Aj Fig. [8] shows that the outage probabilities are quite different in all the five cases as they 
depend on the shape of the region as well as the location of the reference point. We can see that in this 
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Fig. 7. Ratio of the outage probability for tfie reference receiver located at the center and the corner of an L = 3, 4, 5, 6, 7, 8, 9-sided 
polygon for the case of no fading (mo = m = oo) and i.i.d. Rayleigh fading (mo — m = I), with M = 10 interferers i.u.d. in each region, 
reference link distance ro = 5, path-loss coefficient a = 2.5, SINR threshold /3 = dB and SNR po — 20 dB. 

example, the outage probability results for the equilateral triangle, square and disk regions cannot be used 
to provide any meaningful bounds for an arbitrarily- shaped finite region. This re-iterates the importance of 
the proposed frameworks, which allow the outage probability at any arbitrary location of a finite wireless 
network with arbitrary shape to be exactly determined. 



VIII. Conclusion 

In this paper, we have proposed two general frameworks for analytically computing the outage proba- 
bility when a finite number of nodes are distributed at random inside an arbitrarily-shaped finite wireless 
network. For the case that the random nodes are independently and uniformly distributed and the fad- 
ing channels are identically, independently distributed following the Nakagami-m distribution, we have 
demonstrated the use of the frameworks to analytically compute the outage probability at any arbitrary 
location inside the finite wireless network. The probability distribution function of a random node from the 
reference receiver plays a key part in our frameworks and enables the accurate modeling of the boundary 
effects. We have also presented an algorithm to accurately compute the outage probability in a finite 



24 



10 



CO 



g 10 
o 



Q. 


D) 
CC 

■4— ' 

o 



10" 



10" 









^Siit 


Vertex V in Fig. 2 


^^ 


^t 


^ 


\ 


veriex ot iriangie 

Vertex of square 

Circumference of disk 

* Simulations 


\ 


L : 






V-X, ^v 




:^:>^ 


.,_.._ , _■:_■,■_,■_:_,_■,_■,■. 




>S. ', ■" * Hi * 

vv ".v. , T 






..^ ^ 






■^ * -* - -ijt 














] 















-5 







10 



15 



20 



25 



30 



Signal-to-noise ratio, p in dB 



Fig. 8. Outage probability, e, versus tiie signal-to-noise ratio, po, with arbitrary location of reference receiver in the arbitrarily-shaped finite 
region in Fig. 2, with M — 10 interferers i.u.d. in each region, i.i.d. Nakagami-TTi fading channels with mo — m = 3, path-loss coefficient 
a = 2.5, SINR threshold /3 = dB. All the shapes have the same area \A\ = 15000. 



wireless network with an arbitrary shape. A summary of the main outage probability results in this paper 
is shown in Table I 



Appendix A 
Derivation of the Distance Distribution //^(r) 

In this appendix, we summarize the derivation of the distance distribution results using the procedure 



in [[T6|. 

For the uniform BPP, the exact CDF Fji{r) of the distance between a randomly located node and an 
arbitrary reference point located inside a regular convex polygon was derived in p6| . Using this result, 
the PDF of the Euclidean distance between the arbitrary reference point and its zth neighbor node is 



obtained, generalizing the result in [15]. The approach in [16| is also applicable for arbitrarily-shaped 



convex polygons. In this paper, we adapt the procedure in [16| to obtain the PDF //j(r) for an arbitrary 
location of the reference receiver inside an arbitrarily-shaped finite wireless network. The main steps of 
the procedure are summarized below. 
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TABLE III 

Summary main outage probability results. 





Equation 


Node 
distribution 


Region 


Reference receiver 
location 


Fading channels 




m 


any 


any 


any 


any 


•a j^ 
■9 % 


m 


i.i.d. 


any 


any 


i.i.d. 


m 


i.i.d. 


any 


any 


i.i.d. Nakagami-m 


StS 


frg&iu 


i.u.d. 


disk 


any 


i.i.d. Nakagami-m 




([Tg&|3§ 


i.u.d. 


polygon 


center 


i.i.d. Nakagami-m 




e 


any 


any 


any 


Reference link: (|ll|i 


23 


& Interference links: any 


1^ 





i.i.d. 


any 


any 


Reference link: l|ll|( 231 & Interference links: i.i.d. 


P 


in 


i.i.d. 


any 


any 


i.i.d. Nakagami-m with integer mo for the reference link 


^■E 


(|g&{26j 


i.u.d. 


disk 


any 


i.i.d. Nakagami-m with integer mo for the reference link 




@&e 


i.u.d. 


polygon 


center 


i.i.d. Nakagami-m with integer mo for the reference link 



Step 1: Consider a conve)(R polygon with L sides. Let Si and V^ {£ = 1, 2, . . . , L) denote the sides and 
the vertices, which are numbered in an anti-clockwise direction. 

Step 2: Using the property of uniform BPP the CDF Fji{r), which is the probability that the distance 

between an i.u.d. node and the reference receiver is less than or equal to r, is given by 

\ViYo;r)nA\ 0(^0; r) 
^^^'^ = ^4^ = ^AT^ ^''^ 

where V{YQ;r) is a disk region centered at Yq with radius r and 0{YQ;r) is the overlap area between 
V{Yo; r) and the network region, respectively. 



Step 3: In order to find the area of overlap region 0{YQ\r), the approach in [16| is to first find the 



circular segment areas formed outside the sides (denoted as Bi for side t) and the comer overlap areas 
between two circular segments at a vertex (denoted as Ci for vertex t), and then subtract from the area 



of the disk. Thus, ( [35| ) can be expressed as 

FR{r) = 




nr 



J2b,+j2c, 



(36) 



Taking the derivative of ([36]), we have 

fnir) 



^(^-? 



dr 



E 



dc, 

dr 



(37) 



all interior angles are less than tt radians. 
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where d/dr denotes the derivative. Note that -^ and -q^ depend on the radius r and location of Yq. 



Following [16 1, they can be expressed as 



dBf ^ /ps^\ 



2r arccos I — ) , (38) 



dr \ r / 

i —Tc + 6e + arccos ( — - ) + arccos ( — ^^ j j , (39) 



dr 

where ps^ denotes the perpendicular distance from Yq to the side Si, 5e represents the interior angle at 
vertex Ve and the notation Si. Note that for i = I, S^^i = Sl- 

Step 4: In order to find the perpendicular distances ps^, establish a Cartesian coordinate system. For 
arbitrarily-shaped finite regions, without loss of generality, the origin can be placed at vertex Vi. For 
regular convex polygons, the origin can be placed at the center of the polygon to exploit the rotational 



symmetry [16| 



Following the steps above, and substituting the values, the distance distribution /^(r) can be derived. 
This is illustrated in detail in the next subsection. Also the distance distribution results used in Section IVIII 
are summarized. 

A. Reference Receiver Located at Vertex V3 in Fig. |2] 

Consider the arbitrarily- shaped region shown in Fig. |2} with reference receiver located at vertex V3. 

For < r < 7^, the disk region 'D{Yo] r) is limited by sides 5*2 and 53 and vertex V3, i.e., there are two 
circular segment areas outside S2 and S3, which are denoted by B2 and B3. Also there is a comer overlap 
area between them, which is denoted by C3. Hence for this range, /^(r) = r^ (27rr — ^^ — ^^ + ^). 

For IZ < r < \/21Z, the disk region V(Yq] r) is limited by all four sides and three vertices V2, V3 and 



V4, i.e., there are four circular segments formed outside all sides and three comer overlaps between them. 
Hence for this range, /^(r) = r^ f 27rr - f: f^ + E f^V 

For r > \/2TZ, the disk region covers the whole region in Fig. 2 which results in Fji{r) = 1 and 



fnif) = 0. Consequently, in this case, //j(r) is a piece-wise function with two ranges only. Using ( [38| ) 
and ( |39l ), and simplifying, we get ( [32] ) which is used in the generation of the results in Fig. [8j 



27 



B. Reference Receiver Located at Vertex V2 in Fig. |2] 

Using the steps outlined above, the distance distribution fR^r) when the reference receiver is located 
at vertex V2 in Fig. 2 is 

[r. < r < y27^; 



Mr) = ^ 



4 ' 

|r - r arccos (f ) , v^7^ < r < 27^; (40) 




' arccos (— ) — r arccos (— ) , 27^ < t < a/57^; 
The above result is used in the generation of the results in Fig. [8] 

C. Reference Receiver Located at a Vertex in an Equilateral Triangle 

Using the steps outlined above, the distance distribution fR^r) when the reference receiver is located 
at any vertex of an equilateral triangle is 

1 I f r, < r < f7^; 

2r arccos (^) , |7^ < r < y37^; 
The above result is used in the generation of the results in Fig. [8] 

D. Reference Receiver Located at a Vertex in a Square 

Using the steps outlined above, the distance distribution fR^r) when the reference receiver is located 
at any vertex of a square is 

1 f f r, < r < v^7^; 

1^1 fr-2rarccosf ^j , y27^ < r < 27^; 
The above result is used in the generation of the results in Fig. [8] 
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